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Gravitational radiation drives an instability in the r-modes of young rapidly rotating neutron stars. 
This instability is expected to carry away most of the angular momentum of the star by gravitational 
radiation emission, leaving a star rotating at about 100 Hz. In this paper we model in a simple way 
the development of the instability and evolution of the neutron star during the year-long spindown 
phase. This allows us to predict the general features of the resulting gravitational waveform. We 
show that a neutron star formed in the Virgo cluster could be detected by the LIGO and VIRGO 
gravitational wave detectors when they reach their "enhanced" level of sensitivity, with an amplitude 
signal-to-noise ratio that could be as large as about 8 if near-optimal data analysis techniques are 
developed. We also analyze the stochastic background of gravitational waves produced by the r-mode 
radiation from neutron-star formation throughout the universe. Assuming a substantial fraction of 
neutron stars are born with spin frequencies near their maximum values, this stochastic background 
is shown to have an energy density of about 10 -9 of the cosmological closure density, in the range 
20 Hz to 1 kHz. This radiation should be detectable by "advanced" LIGO as well. 

PACS Numbers: 04.30. Db, 04.40.Dg, 97.60.Jd 



I. INTRODUCTION 

Recently Andersson |lj discovered that gravitational 
radiation tends to destabilize the r-modes of rotating 
stars. Friedman and Morsink Q then showed that this 
instability is generic, in the sense that gravitational radi- 
ation tends to make all r-modes in all rotating stars un- 
stable. Lindblom, Owen, and Morsink || have recently 
evaluated the timescales associated with the growth of 
this instability. Gravitational radiation couples to these 
modes through the current multipoles rather than the 
more typical mass multipole moments. This coupling is 
stronger than anyone anticipated for these modes, and is 
so strong in fact that the viscous forces present in hot 
young neutron stars are not sufficient to suppress the 
gravitational radiation driven instability. Gravitational 
radiation is expected therefore to carry away most of the 
angular momentum of hot young neutron stars. These 
results have now been verified by Andersson, Kokkotas, 
and Schutz gj. 

In this paper we study the gravitational waveforms 
that are produced as the r-mode instability grows and 
radiates away the bulk of the angular momentum of a 
hot young rapidly rotating neutron star. The proper- 
ties of the r-modes and the instability associated with 
them are reviewed in Sec. [H]. The equations that de- 
scribe (approximately) the evolution of the r-modes as 
they grow and spin down a rapidly rotating neutron star 
are derived in Sec. These equations are solved nu- 



merically and the results are also presented in Sec. Ill 



The gravitational waveforms associated with the r-mode 
instability are evaluated in Sec. IV. General analytical 



are presented. In Sec. ^ we evaluate the detectability of 
this type of gravitational wave signal by the laser inter- 
ferometergravitational wave detectors such as LIGO || , 
VIRGO ||, and GEO §. We consider the detectability 
of signals produced by single nearby sources, and also the 
detectability of a stochastic background of sources from 
throughout the universe. Finally, in Sec. VI we discuss 



the prospects for gravitational-wave astronomy opened 
up by the r-modes. 



II. THE i?-MODE INSTABILITY 

The r-modes of rotating Newtonian stars are generally 
defined to be solutions of the perturbed fluid equations 
having (Eulerian) velocity perturbations of the form 
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where R and £1 are the radius and angular velocity of the 
unperturbed star, f(r/R) is an arbitrary dimensionless 
function, and Y® m is the magnetic type vector spherical 
harmonic defined by 
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and detailed numerical expressions for these waveforms 



For barotropic stellar models, of primary concern to us 
here, the Euler equation determines the form of these 
modes: The radial dependence f(r/R) is determined 
to be fir / R) — a(r/R) 1 , where a is an arbitrary con- 
stant pj. These modes exist with velocity perturbations 
as given by Eq. ([O]) if and only if I = m ||. Also, the 
frequencies of these modes are given by p] 
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These modes represent large scale oscillating currents 
that move (approximately) along the equipotential sur- 
faces of the rotating star. The restoring force for these 
oscillations is the Coriolis force; hence the frequencies 
of these modes are low compared to the usual / and 
p— modes in slowly rotating stars. These expressions for 
5v and ui are the lowest order terms in an expansion in 
terms of the angular velocity fi. The exact expressions 
contain additional terms of order fl 3 . There may exist 
other modes of rotating barotropic stellar models with 
properties similar to these classical r-modes; however, 
our discussion here is limited to the properties of these 
classical r-modes. 

The density perturbation associated with the r-modes 
can be deduced by evaluating the inner product of v 
(the unperturbed fluid velocity) with the perturbed Eu- 
ler equation, and the equation for the perturbed gravita- 
tional potential Q : 
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The quantity 8^ is proportional to the perturbed gravi- 
tational potential (5$, and is the solution to the ordinary 
differential equation 
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which satisfies appropriate asymptotic boundary condi- 
tions. We note that 8p is proportional to ft 2 and hence 
is small (i.e., higher order in f2) compared to 8v in slowly 
rotating stars. We also note that 8p is proportional to 
YJ_i_i i — having spherical harmonic index one order in I 
highe r th an that of the velocity perturbation. Equa- 
tion (2.4) is the complete expression for the density per- 
turbation to order f2 2 . The exact expression for Sp in- 
cludes additional terms of order il 4 . 

The r-modes evolve with time dependence e 4Wt- */ T as 
a consequence of ordinary hydrodynamics and the influ- 
ence of the various dissipative processes. The real part 
of the frequency of these modes, u>, is given in Eq. ((2.3[), 
while the imaginary part 1 /r is determined by the effects 
of gravitational radiation, viscosity, etc. The simplest 
way to evaluate 1 jr is to compute the time derivative of 
the energy E of the mode (as measured in the rotating 
frame) . E can be expressed as a real quadratic functional 
of the fluid perturbations: 



E = l 
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Thus the time derivative of E is related to the imaginary 
part of the frequency 1 /r by 
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Since the specific expressions for the time derivative of E 
due to the influences of gravitation al r adiation 1 1 and 
viscosity [|l2| are well known, Eq. (2.7) may be used to 
evaluate the imaginary part of the frequency. 
It is convenient to decompose 1/t: 
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where 1/t gr , 1/t s , and l/r B are the contributions due to 
gravitational radiation emission, shear viscosity and bulk 
viscosity respectively. Expressions for these individual 
contributions for the r— modes are given by [pf: 
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and 

1 AR 21 - 2 
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where Sp is given in Eq. (2.4). We note that the ex- 
pression for 1/t b in Eq. (2.11) is only approximate. The 



exact expression should contain the Lagrangian density 
perturbation Ap in place of the Eulerian perturbation dp. 
The bulk viscosity (see Eq. [2.15 ) is a very strong function 
of the temperature, being proportional to T 6 . Thus, the 
result of any error that might occur in our approxima- 
tion for l/r B is simply to shift slightly the temperature 
needed to achieve a given viscosity timescale. Numeri- 
cal estimates show that changing this quantity by even 
a factor of one hundred (as suggested by Ref. ) , docs 
not substantially affect the important physical quantities 
computed here (i.e. the spindown rate or the final angular 
velocity of the star). 

We have evaluated these expressions for the imaginary 
parts of the frequency for a "typical" neutron star model 
with a polytropic equation of state: p = kp 2 , with k cho- 
sen so that a 1.4M model has the radius 12.53km. We 
use the usual expressions for the viscosity of hot neutron 
star matter 13 1: 
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where all quantities are expresse d in c gs un its. We have 
evaluated the expressions, Eqs. ( |2.9| ) (2.11), for the dis- 
sipative timescales with fiducial values of the angular ve- 
locity £1 = \JixGp and temperature T = 10 9 K. These 
fiducial timescales f GR , f v , and f B are given in Table || 
for the r-modes with 2 < I < 6. It will be useful in 
the following to define a timescale associated with the 
viscous dissipation l/r v = 1/t s + l/r B . The viscous 
timescale t v and the gravitational timescale 

7~gr can be 

expressed then in terms of the fiducial timescales in a 
way that makes their temperature and angular velocity 
dependences explicit: 
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where / is the moment of inertia of the equilibrium state 
of the star, and J c is the canonical angular momentum 
of the r-mode. 

In this simple model of the star the angular momen- 
tum J is a function of the two parameters that charac- 
terize the state of the system: J = J(fi,a). We can 
determine this functional relationship approximately as 
follows. The canonical angular momentum of an r-mode 
can be expressed in terms of the velocity perturbation Sv 

byR, 



Jr. 



I 



2(u + m) 



pSv ■ 5v*d x 



(3.2) 



For the I = 2 r-mode of primary interest to us here this 
expression reduces (at lowest order in £1) to 



J c = -ma 2 J MR 2 , 



(3.3) 
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where J is defined by 
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III. EVOLUTION OF THE 7?-MODES 

To determine the gravitational waveform that will re- 
sult from the instability in the r-modes, we must estimate 
how the neutron star evolves as the instability grows and 
radiates the angular momentum of the star away to in- 
finity. Initially the mode will be a small perturbation 
that is described adequately by the linear analysis that 
we have described above. However, as the mode grows, 
non-linear hydrodynamic effects become important and 
eventually dominate the dynamics. At the present time 
we do not have available the tools to follow exactly this 
non-linear phase of the evolution. Instead, we propose a 
simple approximation that includes (we believe) the basic 
features of the exact evolution. 

We treat the star as a simple system having only two 
degrees of freedom: the uniformly rotating equilibrium 
state parameterized by its angular velocity £7, and the 
r-mode parameterized by its amplitude a. The total an- 
gular momentum J of this simple model of the star is 
given by, 



j = m + J c 



(3.1) 



TABLE I. Gravitational radiation and viscous timescales 
(in seconds) are presented for T = 10 9 K and Q, = ^/nGp. 
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TGR 



T S 



Tb 



-3.26 x 10 u 
-3.11 x 10 1 
-2.85 x 10 2 
-2.37 x 10 3 
-1.82 x 10 4 



2.52 x 10 s 
1.44 x 10 s 
1.07 x 10 s 
8.79 x 10 7 
7.58 x 10 7 



6.99 x 10 s 
5.13 x 10 s 
4.01 x 10 s 
3.26 x 10 s 
2.74 x 10 s 



For the polytropic models studied in detail here the di- 
mensionless constant J = 1.635 x 10~ 2 . The moment of 
inertia I can also be conveniently expressed as 
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where / is given by 



I = 



8tt 



MIR 2 J 



pr dr. 



(3.5) 



(3.6) 



For the polytropic models considered here I — 0.261. 
Thus, our simple model of the angular momentum of the 
perturbed star is 



J(n,a) = (I- iJa 2 )flMR 2 . 



(3.7) 



The perturbed star loses angular momentum primar- 
ily through the emission of gravitational radiation. Thus, 
we compute the evolution of J(Q, a) by using the stan- 
dard multipole expression for angular momentum loss. 
The I = 2 r-mode is the primary source of gravitational 
radiation in our simple model of this system, and this 
mode loses angular momentum primarily through the 
I = m = 2 current multipole. Thus the angular mo- 
mentum of the star evolves as 
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The I = m = 2 current multipole moment S22 for this 
system is given by 
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Combining Eq. ( |3.8| ) for t he angula r mo mentum evolu- 
tion of the star with Eqs. (^), and (3/7), we obtain one 
equation for the evolution of the parameters fl and a that 
determine the state of the star: 
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During the early part of the evolution of the star, the 
perturbation analysis of the r-modes described earlier ap- 
plies. In addition to radiating angular momentum from 
the star via gravitational radiation, the mode will also 
lose energy via gravitational radiation and neutrino emis- 
sion (from the bulk viscosity) and also deposit energy into 
the thermal state of the star due to shear viscosity. It is 
most convenient to obtain the equation for the energy 
balance during this part of the evolution in terms of the 
energy E of the mode as defined in Eq. (2.6). For the 
I = 2 r-mode E is given by 
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The time derivative of E is precisely the quantity that 
was used to determine the im aginary part of the fre- 
quency of the mode in Eq. ( |2.7| ): 
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Equation (3.12) together with ( |3 . 1 1 ) therefore provides 
a second equation for determining the evolution of the 
parameters 51 and a that specify the state of the star: 



da dn / 1 

fi— + a— = -afl 

dt dt \ t gr 
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Equations (3.1C) and (3.13) can be combined then to 
determine the evolution of Q and a during the portion of 
the evolution in which the perturbation remains small: 
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The equation of st ate de pend ent parameter Q that ap- 
pears in Eqs. fl3.14| ) and ( j3J5| ) is defined by Q = 3 J/21. 
For the polytropic model considered in detail here Q = 
9.40 x 10~ 2 . We note that during the initial linear evo- 
lution phase the angular velocity of the s tar Q is nearly 
constant, evolving according to Eq. ( 3.14| ) on the viscous 
dissipation timescale. During this phase the amplitude 
of the mode a grows exponentially on a timescale that is 
comparable to the gravitational radiation timescale. 

After a short time (about 500 s in our numerical so- 
lutions) the amplitude becomes so large that non-linear 
effects can no longer be ignored. We have not yet de- 
veloped the tools needed to follow the evolution exactly 



during this non-linear phase. However, we do have some 
intuition about the non-linear hydrodynamical evolution 
of gravitationally driven instabilities in rotating stars. 
This intuition comes from the studies of the effects of 
gravitational radiation reaction on the evolution of the 
ellipsoidal models Jl5|Jl6| . In that case the unstable mode 
grows exponentially until its amplitude is of order unity. 
At that point a kind of non-linear saturation occurs, 
and the growth of the mode stops. The excess angu- 
lar momentum of the star is radiated away and the star 
evolves toward a new lower angular momentum equilib- 
rium state. We expect a similar situation to pertain in 
the evolution of the r-modes. Thus, we expect non-linear 
effects will saturate and halt the further growth of the 
mode when the amplitude of the mode becomes of order 
unity. Thus, when the amplitude a grows to the value 

a 2 = k, (3.16) 

(where k is a consta nt of order unity ) we stop evolving 
the star using Eqs. (3.14) and (3. IE). Instead we set 
da I ' dt = during the saturated non- linear phase of the 
evolution, while continuing to evolve the angular velocity 
f2 by Eq. ( 3.10 ) as angular momentum is radiated away 
to infinity by gravitational radiation. During this phase, 
then the angular velocity evolves by 



d£l _ 2Q kQ 
dt t G r 1 - kQ 



(3.17) 



The r-mode will evolve during the saturated non- 
linear phas e of it s evo lution approximately according to 
Eqs. ( 3.16 ) and ( 3.1 7| ) . During this phase the star will 
lose most of its angular momentum, and spin down to 
a state having an angular velocity that is much smaller 
than Qk ss Zy/nGp. The star will eventually (in about 1 
year in our numerical solutions) evolve to a point where 
the angular velocity and temperature become sufficiently 
low that the r-mode is no longer unstable. The end of the 
evolution is characterized by a phase in which the viscous 
forces and gravitational radiation damp out the energy 
remaining in the mode and move the star slowly to its fi- 
nal equilibrium configuration. During this final phase, 
the mode is again of small amplitude and so the lin- 
ear approximation is adequate to describe the evol ution . 
We monitor the quantity on the right side of Eq. (3.15) 
throughout the non-linear evolution phase. When it be- 
comes nega tive w e cha nge t he evolution equations again, 
from Eqs. ( 3.16 ) and ( 3.17 ) back to the linear equations 
Eqs. (fr|) and ( |3~T5| ). 

In summary then, we model the evolution of the r- 
mode as having three distinct phases: (i) The hot young 
neutron star is born rapidly rotating with a small ini- 
tial excitation in the I = 2 r-mode. T his m ode in itiall y 
grows exponentially according to Eqs. ( 3.14| ) and ( 3.15 ). 
(ii) The amplitude of the mode saturates due to non- 
linear hydrodynamic effects at a value of order unity. 
The bulk of the angular momentum of the star is ra- 
diated away by gravitational radiation during this phase 
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FIG. 1. Angular velocity evolution as gravitational radia- 
tion spins down a hot rapidly rotating neutron star, a mea- 
sures the amplitude of the initial perturbation. The dashed 
curve shows the critical angular velocity Q. c above which the 
r-modes are unstable. 



according to Eqs. ( 3.16 ) and ( 3.17 ). (iii) The final phase 
of the evolution begins when the right side of Eq. ( 3.15| ) 
becomes negative so that the mode begins to be damped 
out. During the f inal p hase t he star evolves again accord- 
ing to Eqs. ( fig ) and (flj|). 

In order to complete our model for the evolution of 
the r-modes we must specify how the temperature of the 
star evolves with time. We do this by adopting one of the 
standard descriptions of the cooling of hot young neutron 
stars. These stars are expected to cool primarily due to 
the emission of neutrinos via a modified URCA process. 
The temperature during this phase falls quickly by a sim- 
ple power law cooling formula §L7\ : 



T(t) 
10 9 K 



I ( 10 9 K 
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Ti 



-1/6 



(3.18) 



where Ti is the initial temperature of the neutron star, 
and t c is a parameter that characterizes the cooling rate. 
For the modified URCA process r c w ly. A typical value 
for the initial temperature is Tj « 10 11 K . Equ ation ( 3. IS ) 
can now be inserted into Eqs. (3.14)-( |3.17 ) to provide 
explicit differential equations for the time evolution of 
the angular velocity of the star and the amplitude of the 
mode. These equations can be solved numerically in a 
straightforward manner. 

Figs, [j] and || illustrate the solutions to these equations. 
The dashed curves in Figs. |l| and || show the critical an- 
gular velocity fl c , defined by l/r(O c ) = 0, above which 
the r-modes are unstable. Figure |l| shows the evolution 
(f2 plotted versus T) of the angular velocity of the star 
for k = 1.0 and a range of values of the initial value of 
the parameter a. In these simulations we have assumed 
that the initial angular velocity of the star is f2 = VIk- 
This figure illustrates that the final non- linear part of the 
evolution is remarkably insensitive to the initial size of 
the perturbation. 

Figure || illustrates the dependence of the evolution on 
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FIG. 2. Evolution of the angular velocity of the star de- 
pends on the parameter k in the non-linear saturated phase, 
but the final angular velocity of the star is insensitive to this. 
The dashed curve is the same as in Fig. 1. 



the parameter k by showing several evolutions with ini- 
tial values of a = 10~ 6 . The parameter k measures the 
degree of saturation that occurs in the non-linear spin- 
down phase of the star. In our numerical studies we ex- 
amine the limited range 0.25 < k < 2. If n is taken to be 
too small, the mode simply does not grow to the point 
that non-linear effects can stop its growth. Conversely 
if k is taken too large, then our simple evolution equa- 
tions based in part on t he lin ear perturbation theory be- 
come singular, e.g. Eq. ( 3.17 ). The equations break down 
because the (negative) canonical angular momentum of 
the mode equals the (positive) angular mome ntum of the 
equilibrium configuration, and therefore Eq. (3.1) yields 
the unphysical result that the star has no net angular 
momentum. 

We artificially stop all of our evolution curves when the 
temperature of the star falls to 10 9 K. Below this tem- 
perature we expect superfluidity and perhaps other non- 
perfect fluid effects to make our simple simulation highly 
inaccurate |l8j . Figure || illustrates that the gravitational 
radiation instability in the r-modes is nevertheless effec- 
tive in radiating away most of the angular momentum of 
the star before the star cools to the point that superflu- 
idity or other effects are expected to become important. 
Figure |^ shows that the amount of angular momentum 
lost in this process is remarkably insensitive to the value 
of k. Thus, the final upper limit on the angular velocity of 
the star is (fortunately) fairly insensitive to our assump- 
tion about the exact nature of the non-linear portion of 
the star's evolution. 

In this simple model of the evolution of the unsta- 
ble star, we have ignored the effect that viscous heating 
might have on the cooling rate of the star. If there were 
too muc h vis cous heating, then the cooling formula given 
in Eq. (3.18) would not be correct. We have evaluated 
the importance of this re-heating effect by comparing the 
rate at which thermal energy is being radia ted away from 
the star by neutrinos according to Eq. ( 3. 18| ) with the rate 
that viscous dissipation deposits thermal energy into the 
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FIG. 3. Cooling rate due to neutrino emission is compared 
to the viscous heating rate in our numerical evolution of the 
unstable r-modes. 
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FIG. 4. Angular velocity evolution due to gravitational ra- 
diation emission using several values for the timescale r c on 
which the star cools. 



star. Neutrino cooling removes energy from the thermal 
state of the star at the rate [Il7) 



dU 
~dt 



7.4 x 10 39 ( 
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ergs/s. 



(3.19) 



Thermal energy is generated by shear viscosity as the star 
evolves, but not by bulk viscosity. Bulk viscosity radiates 
away its excess energy directly by neutrinos without sig- 
nificantly interacting with the thermal energy contained 
in the star. Thus, energy is transfered from the canonical 
energy of the r-mode to the thermal energy of the star 
by the formula 



dE c _ 2a 2 £l 2 MR 2 J 
dt t s 
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Figure |^ compares the values of dU/dt and dE c /dt 
for our numerical evolution with k = 1.0 and initial 
a = 10~ 6 . We find that viscosity does not deposit en- 
ergy into the thermal state of the star at a significant 
rate until the temperature of the system falls to about 
10 9 K. At this point the star has already lost most of its 
angular momentum to gravitational radiation, and other 
dissipative effects (such as those associated with super- 
fluidity) which are not modeled here will start to play 
a significant role jlS!. Thus, we are justified in ignoring 
the effects of viscous re-heating on the thermal evolution 
of the star during the early part of its evolution modeled 
here. 

The modified URCA process that determines the ther- 
mal evolution used in our evolutions is the standard 
mechanism by which neutron stars are expected to cool 
down to about 10 9 K. Other less standard mechanisms 
have also been proposed which could significantly speed 
up the cooling pj^ . These mechanisms include neutrino 
emission processes that require the presence of exotic 
species (such as quarks or pions) as free particles in the 
cores of these stars. We have ignored these possibilities 
in the evolutions described above. If these particles do 
exist in the cores of neutron stars, we expect that they 



will only be present in a small volume of material at the 
centers of these stars. The cores of these stars may well 
cool rapidly, but the outer layers where the r-m ode i s 
large will continue to cool at the rate given in Eq. ( 3.1§| ) 
until thermal conduction can move energy from the outer 
layers back into the core. 

To estimate what effect a somewhat more rapid cool- 
ing might have on the evolution of the r-modes, we have 
artificially varied the value of the parameter r c that deter- 
mines the cooling rate in Eq. (3.18). Figure ^ shows the 
results of the evolution of an r-mode with initial a = 10~ 6 
and k = 1.0 for several values of r c . We see that while the 
details of the evolution are somewhat effected, the total 
amount of angular momentum radiated away by gravita- 
tional radiation, and the final angular velocity of the star 
are fairly insensitive to the rate at which the star cools 
to 10 9 K. 



IV. GRAVITATIONAL WAVEFORMS 

As the r-mode grows and evolves it emits gravitational 
radiation. In this section we calculate the waveforms for 
the gravitational wave strain h{t) and its Fourier trans- 
form h(f) that are produced by this r-mode instabil- 
ity. These are the quantities that can be measured by 
the gravitational wave detectors now under construction 
(LIGO, VIRGO, GEO, etc.). During the non-linear sat- 
uration phase of the r-mode evolution, gravitational ra- 
diation controls the dynamics. In this case h(f) turns 
out to be independent of the details of the evolution, and 
thus can be determined quite generally. An evolutionary 
model is needed, however, to determine df /dt, the time 
dependence of the various quantities, and the initial and 
final frequencies of the saturation phase. In this section 
we present the general model independent derivation of 
h(f) for the non-linear saturated phase of the evolution. 
We also give expressions for the gravitational wave strain 
that apply to the early phases of the evolution using the 
simple model discussed in Sec. III. 
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The frequency-domain gravitational waveform 

h(f)= / e 2mft h(t)dt (4.1) 



is determined completely by the assumption that the an- 
gular momentum radiated as gravitational waves comes 
directly from the angular momentum of the star. This as- 
sumption is expected to be satisfied during the non- linear 
saturated phase of the evolution, but not during the early 
evolution when the mode is growing exponentially. This 
derivation is based on Blandford's analysis (as discussed 
in jnj) of white dwarf collapse to a neutron star which is 
halted by centrifugal forces (see also [^,^l) ) . Such a star 
can only collapse to a neutron star by shedding its ex- 
cess angular momentum through gravitational waves. In 
that situation as in the non-linear saturation phase of the 
r-mode evolution, gravitational radiation determines the 
rate at which angular momentum leaves the system, and 
this in turn determines the rate at which the frequency 
of the radiation evolves with time. 

In the stationary phase approximation (which is al- 
ways valid for a secular instability) the gravitational wave 
strain h(t) is related to its Fourier transform h(f) by 



\h(t)\ 2 = \h(f)f 



(4.2) 



Throughout this discussion we treat h{t) as a complex 
quantity with purely positive frequency. For the I = 2 
mode of primary importance here, the mode frequency 
is uj = or f — where / is the frequency of the 

emitted gravitational waves measured in Hz. Assuming 
the star is uniformly rotating, its angular momentum is 
J = in, where I is the star's moment of inertia. The 
moment of inertia is fairly independent of angular veloc- 
ity (especially at small angular velocities where most of 
the detectable signal from these sources is likely to orig- 
inate) and is also fairly independent of the amplitude of 
the excited r-mode (see Eq. 3.7). Thus, I is reasonably 
well approximated by its non-rotating value. Thus 



dj 37T 



(4.3) 



The rate at which angular momentum is radiated away 
by a source is related to the gravitational wave amplitude 
by the expression [O 



dJ ~m 1 2 2 

— = 4tt.D u (hi. 

dt uj 16tt x + 



(4.4) 



where h + and h x are the amplitudes of the two polar- 
izations of gravitational waves, D is the distance to the 
source, and (. . .) denotes an average over the orientation 
of the source and its location on the observers sky. U sing 
dt/df = dJ/dfidJ/dt)- 1 and combining Eqs. (gg)-(|j) 
we obtain 



M/)I 2 + IM/)I 2 > = 



3/ 



(4.5) 



The measured value of \h(f)\ 2 depends on the orienta- 
tion of the source and its location on the detector's sky. 
Averaged over these angles, its value is given by 



(IM/)I 2 ) = *(IM/)I 2 + IM/)I 2 



(4.6) 



We are actually interested in the average over source loca- 
tions in three-dimensional space, not just the two angles 
on the sky. The spatial average weights more strongly 
those orientations that yield stronger signals, effectively 
increasing (\h(f)\ 2 ) by about |. Combining these re- 
sults then, the average value of h produced by our fidu- 
cial r-mode source (with M — 1.4M Q , D — 20Mpc, 
R = 12.5km) is 



h(f) = 5.7 x 10 



-25 



'1 kHz 



/ 



Hz~ 



(4.7) 



Note that this expression does not depend (in the fre- 
quency domain) on the details of the evolutionary model 
apart from the u ppe r and lower frequency cutoffs. This 
expression, Eq. ( |4.7| ), only depends on the assumption 
that the frequency of the mode evolves as angular mo- 
mentum is radiated by the star according to Eq. ( [1.3] ). 
We expect this to be satisfied during the non-linear sat- 
urated phase of the r-mode evolution, but probably not 
during the early linear growth phase of the mode. 

To obtain the complete waveforms for a particular evo- 
lutionary model, we start with the usual expression for 
the gravitational field in terms of its multipoles |ll[] . We 
average this expression over angles in the manner de- 
scribed above to obtain 



hit) 



3 uj 2 S 2 
80tt D 



(4.8) 



For the simple two-parameter evolution of the r-mode 
instability described in S ec. Ill, this expression can be 
simplified using Eq. (3.9) to 



h(t) = 4.4 x 10~ 24 a 



20Mpc 
D 



(4.9) 



This simple evolutionary model also gives a simple for- 
mula for the frequency evolution. Dur ing t he non-linear 
saturation phase of the evolution Eq. (3.17) can be writ- 
ten as 



df_ 
dt 



-1.8k 



1 kHz 



Hz/s, 



(4.10) 



where we have assumed that kQ <C 1. The time for the 
gravitational wave frequency to evolve to its minimum 
value / m in is obtained by integrating Eq. (4.10): 



t 



1.0 /l20Hz\ 

K \ /min / 



(4.11) 



Analogous model dependent expressions can also be 
derived for the early linear phase of the evolution. During 
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this period the amplitude a grows expo nenti ally on the 
gravitational timescale according to Eq. ( 3.14 ), while the 
frequency of the mode changes extremely slowly accord- 
ing to Eq. (3.15). Solving these equations approximately 



gives 



4f 
dt 



2.7a 2 



f 



1 kHz 



(4.12) 



Using Eq. (4.9) this implies that h during the linear 
growth phase is given approximately by 



h{f) w 4.7 x lQ- 25 Vt 



f 



1 kHz 



3/2 



Hz" 



(4.13) 



where t = t(f) is obtained from Eq. (4.12). The factor 
pl 2 that appears in Eq. (4.13) is essentially constant, 



being given by the initial mode frequency as determined 
by the initial angular velocity of the star. Since this 
factor is essentially constant during the linear evolution 
phase, this implies that h(f) grows as y/i. The duration 
of the initial linear growth phase can be estimated from 
the solution for the amplitude a: 



At = 1.5 x 1(T 



1 kHzV 



In Gv/k/oq) 
ln(10 6 ) 



(4.14) 



where a is the initial size of the perturbation. Dur- 
ing this interval, the mode frequency decreases by an 
amount of order 0.1k Hz. This is the width of the initial 
spikes shown in Figs. 7 and 8. The maximum amp litude 
achie ved b y h can also be determined from Eqs. ( 4.12 ) 
and (fllj): 



max(/i) 1.8 x 10 



-23 



)/2 



1 kHz\ 

log (■</*,/ a ) 



log(10 6 ) 



1/2 



Hz" 



(4.15) 



This expression for max(/i) is fairly insensitive to the du- 
ration of the growth phase, and as well to the exact point 
at which the transition to the non-linear saturation phase 
occurs. The value of max(/i) is rather sensitive to our ex- 
pression for df /dt however. If df /dt were to differ during 
the linear growth phase from the expression used here by 
any small effect (such as non-linear modifications of the 
frequency of the mode) then the the resulting change on 
max(/i) co uld b e l arge. The analytical expressions given 
here, Eqs. (4.13)— (4.15), do however accurately represent 
(to within a few percent) the exact numerical solutions 
to the equations for our simple model in this region. 

Figure |E| illustrates our full numerical solutions for the 
time dependence of the gravitational wave amplitude h(t) 
for several values of the parameters a with k = 1 . Fig- 
ure H illustrates the time dependence of h(t) for various 
values of n with a = 10~ 6 . All of these curves represent 
the gravitational radiation emitted by a neutron star ini- 
tially spinning with angular velocity f2 = Qk — \\I^Gp. 
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FIG. 5. Time dependence of the gravitational wave ampli- 
tude h(t) for a detector located at D — 20Mpc. The peak 
amplitude is very insensitive to the initial size of the pertur- 
bation Q. 
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FIG. 6. Time dependence of the gravitational wave ampli- 
tude h(t). The late time amplitude is rather insensitive to the 
non-linearity parameter k. 



In this section and the next we terminate the evolution 
once the star has cooled to 10 9 K. Below this temperature 
the evolution will be significantly affected by mechanisms 
not dealt with in this paper, such as the superfiuid tran- 
sition, the re-heating of the star by viscous dissipation 
in the mode, and dissipation mechanisms (e.g. plate tec- 
tonics) associated with the rapidly forming crust. 

Figs. and || illustrate the frequency dependence of 
h(f) for a 1.4M Q neutron star located at 20Mpc using 
the same values of a and k used in Figs. || and ||. Fig- 
ure [?] illustrates that h(f) is remarkably insensitive to 
the initial size of the perturbation a. The sharp verti- 
cal spikes appearing at the high-frequency ends of these 
curves are due to the extremely monochromatic gravi- 
tational waves emitted during the linear growth phase. 
The structure of this sp ike in our model is accurately de- 
scribed by Eqs. ( [1.13| )-( 4.15 ), but it is not clear that this 
spike is a robust feature of our model. During the phase 
of the evolution that produces the spike, the amplitude 
of the mode is quite small except for a period of about 
one minute. Thus the total amount of radiated energy 
and angular momentum contained in this spike is quite 
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FIG. 7. Frequency dependence of the gravitational wave 
amplitude h for a source located at 20 Mpc. 
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FIG. 8. Frequency dependence of the gravitational wave 
amplitude h for various values of the initial angular velocity 
of the star. 

small. The spike is not likely to play an important role in 
the detection of these sources even if it is a real feature 
of the r-mode instability. 



V. DETECTABILITY 

In this section we discuss the detectability of the grav- 
itational radiation emitted by young neutron stars spin- 
ning down due to the r-mode instability. This radiation 
might be detected as strong sources from single spindown 
events or as a stochastic background made up of many 
weaker sources. 



A. Single Sources 

First we estimate the signal-to-noise ratio S/N for a 
single source located at a distance D, chosen to be large 
enough so that there is likely to be a reasonable event 
rate (say a few events per year). This distance D must 
be large enough then to include several thousand galaxies 
(assuming that the observed supernova rate is compara- 
ble to the neutron star formation rate). Thus we take 



this fiducial distance to have the value D = 20 Mpc p2| , 
the approximate distance to the Virgo cluster of galaxies. 
We estimate the optimal value of S/N that could be ob- 
tained by matched filtering. Because matched filtering is 
probably not feasible for these sources, this estimate pro- 
vides only an upper limit of what might be achieved. We 
briefly discuss two more realistic search strategies based 
on barycentered Fourier transforms of the data ]2q| . 

Using matched filtering, the power signal-to-noise ratio 
(S/N) 2 of the detection is given by 



IKfW 

Sh(f) 



df 



(5.1) 



where Sh (/) is the power spectral density of the detector 
strai n n oise. The constant in front of the integral in 
Eq. fl4.5p is 2 (instead of 4 as in Ref. p4|) because our h is 



complex (with purely positive frequency). Equation (5.1) 
can also be written 



S 



'1 
f 



(5.2) 



N J 'J f \hr* 

Here the rms strain noise h lms in the detector is given by 



VfShU), 



(5.3) 



where Sh(f) is the power spectral density of the noise, 
and the characteristic amplitude h c of the signal is de- 
fined by pi 



(5.4) 




The quantity multiplying h on the right side of Eq. (5.4) 
is generally interpreted as the number of cycles radiated 
while the frequency changes by an amount of order /. 
This interpretation is correct as long as the frequency 
evolution is very smooth. 

However, our evolutionary model of the frequency evo- 
lution contains a discontinuity as the mode stops linearly 
evolving and saturates. Consequently the actual num- 
ber of cycles spent near the initial frequency is far fewer 
than indicated by Eq. (5.4). The quantity h c is a useful 
estimator of the effective filtered amplitude of a signal 
because the integral of (h c /h lms ) 2 always gives the opti- 
mal (S/N) 2 . Therefore a spike in h c must be interpreted 
with some caution — the peak value of the spike is useless 
unless one also knows the bandwidth of the spike. 

In Figure || we plot h c versus frequency, superimposed 
on h Tms for three LIGO interferometer configurations. In 
the saturation phase (i.e., not including the spike) h c is 
well approximated by 



5.7 x 10 
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(5.5) 
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FIG. 9. Characteristic gravitational wave amplitude h c 
(solid curve) compared to the noise amplitude /i rms (dashed 
curves) for LIGO. 



We plot h lms for the LIGO "first interferometers" || 
(which we abbreviate LIGO I), the "enhanced interfer- 
ometers" pq] (LIGO II), and the "advanced interferom- 
eters" H (LIGO III). The noise power spectral density 
for LIGO I is well approximated by the analytical fit |2(| 



S h (f) = y 



f 



fo 



(5.6) 



where S Q = 4.4 x 10~ 46 Hz" 1 and f Q = 175 Hz. For 
LIGO II we construct the approximation 



9/2 



1+ 4- 



(5.7) 



where S a = 8.0 x 10~ 48 Hz" 1 and f Q = 112Hz. For 
LIGO III the noise spectral density is well approximated 

by p7| 



Sh(f) = 




(5.8) 



where S D = 2.3 x 10~ 48 Hz" 1 and f Q = 76 Hz. 

Most of the contribution to S/N in Eq. (5.2) comes 
from the saturation phase of the evolution which is 
largely model independent as discussed earlier. Given 
a detector noise curve, S/N is thus independent of most 
details of the waveform except the final frequency of the 
neutron star. We therefore have 



9/ 
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(5.9) 



The minimum frequency / m i n reached by the r-mode evo- 
lution is about 120 Hz. At frequencies slightly larger than 
/ m in, the LIGO II noise is dominated by photon shot 
noise. If we ignore the other noise components, Sh(f) 
becomes 



S h (f) = 2.6 x 10" 52 / 2 



(5.10) 



For / max ^> /min the integral is dominated by the lower 
cutoff. Thus, the LIGO II S/N is given approximately 
by 



S_ 

N 



— « 8.8 
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10 45 cgs 



1/2 
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fn 



120 Hz 



(5.11) 



Including the other components of the noise decreases 
S/N somewhat. For the numerical evolutions of our sim- 
ple model with n = 1.0 which terminate at 10 9 K, we 
find we find S/N = 1.2, 7.6, and 10.6 for LIGO I, II, and 
III respectively at D = 20 Mpc. These last three results 



scale with / and D just as in Eq. (5.11), but the depen- 
dence on /min is now more complicated. The contribution 
to S/N from the high frequency spike in our model is 0.6 
for LIGO II, and about 0.1 for LIGO I. While the height 
of the spike in h c may not be a robust feature of our sim- 
ple model, the contribution that this spike makes to the 
overall S/N probably is. These numbers indicate that 
the gravitational radiation from the r-mode instability is 
a probable source for LIGO II if some near-optimal data 
analysis strategy can be developed. It appears unlikely 
that the radiation from the high frequency spike will be 
detectable. 

Given the recent discovery of the ultrafast young pul- 
sar PSR J0537-6910 in the Large Magellanic Cloud ||, 
it is interesting to examine the effect on dctectability of 
a relatively high superfluid transition temperature. If 
the initial period of PSR J0537-6910 was about 7 ms (as 
extrapolated from the braking indices of typical young 
pulsars) , that could indicate a superfluid transition tem- 
perature of about 2 x 10 9 K and final gravitational-wave 
frequency of about 200 Hz. Cutting off our simulations 
at this point in the evolution, we obtain for LIGO II the 
S/N of about 5. 

It is clear from Fig. || that the first-generation LIGO 
and VIRGO detectors will not see r-mode events from 
the Virgo cluster. Their sensitivity is a factor of about 
8 worse than the enhanced detectors. We have also con- 
sidered the possibility that GEO might detect these 
sources by using narrow-banding, where it can improve 
its sensitivity in a restricted frequency range at the ex- 
pense of worse sensitivity elsewhere. However, for the 
kind of broad-spectrum signal produced by the r-mode 
instability, narrow-banding is in fact neutral: the gain of 
signal-to-noise ratio in the selected band is just compen- 
sated by the loss over the rest of the spectrum. So GEO 
is not likely to see these signals either. Nor will advanced 
resonant detectors: at their frequencies and in their rela- 
tively narrow bandwidths, there is just not enough power 
in these signals if the sources are in the Virgo cluster. For 
example, the proposed GRAIL detector Q operating at 
its quantum-limited sensitivity (Sh = 1.6 x 10 -48 Hz -1 ) 
between 500 and 700 Hz would have S/N w 1 for a source 
at the distance of the Virgo cluster. 

Matched filtering using the year-long waveform tem- 
plates that would be needed to describe these sources 
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completely could yield the signal-to-noise ratios quoted 
above. However, this is not a practical strategy for this 
type of signal due to the prohibitively large number of 
templates that would be needed to parameterize our ig- 
norance of these sources and the resulting high compu- 
tational cost of filtering the data with these templates. 
Other strategies equivalent to combining the results of 
shorter template searches might well be computationally 
feasible, although they would obtain less than the opti- 
mal S/N. 

The barycentered fast Fourier transform (FFT) tech- 
nique that has been designed to search for nearly pe- 
riodic signals ]23j] might well provide one such method. 
Figure 10 shows the spindown age 



Tsd 



~ ! df 



(5.12) 



for a neutron star spinning down due to the r-mode insta- 
bility. This quantity provides a reasonably good estimate 
of the time spent by the evolving star in the saturated 
non-linear phase, but it is not a good estimate of the 
amount of time spent during the linear growth phase (for 
the reasons outlined above). During the saturated phase 
T s( j is given approximately by 
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The signal becomes quite monochromatic after about one 
day, in the sense that the spindown timescale is long com- 
pared to the inevitable daily modulation of the signal due 
to the motion of the Earth. Thus the search techniques 
for periodic sources should work well. In many cases the 
supernova will be observed, yielding the location of the 
source and allowing a search over spindown parameters 
only, a directed spindown search. 

The most straightforward way to conduct a directed 
spindown search is to search over generic spindown pa- 
rameters as discussed by Brady and Creighton j3(| . This 
involves re-sampling the data so as to render sinusoidal a 
signal with arbitrary (smooth) frequency evolution, then 
Fourier transforming it. The signal frequency evolution 
(which determines the re-sampling) is modeled as 



N 1 

/(t) = /o£- 



t 



(5.14) 



where /o is the frequency at the beginning of the FFT, 
t is the time measured from the beginning of the FFT, 
and T n are expansion parameters with t\ — T s d- The 
number N of spindown parameters needed is set by the 
requirement that the frequency drift due to the next term 



in the series (5.14) be less than one frequency bin of the 



FFT (which is in turn determined by the integration time 
Tint). This implies (assuming r n ~ T s d) that 
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One must choose a set of points in the spindown 
parameter space for which to perform the re-sampling 
and FFT. Too few points and one misses signals by re- 
sampling at the wrong rate; too many and the compu- 
tational cost of performing all the FFTs becomes pro- 
hibitive if data analysis is to be performed "on-line," i.e. 
keep up with data acquisition. The points in parameter 
space are chosen using a metric which relates distance 
in parameter space to loss of S/N p3[ |. The integration 
time, which is set so as to optimize the sensitivity of an 
on-line search for the (fixed) computational power avail- 
able, is far shorter than a year even for a teraflop com- 
puter. Therefore such a search would achieve at best a 
fraction of the optimal S/N . It is possible that a hi- 
erarchical version of this strategy could be developed, in 
which the best candidates from a year of shorter FFTs 
are somehow combined to give an improved confidence 
level. Developing such a strategy would require exten- 
sive further investigation. 

One way to increase the sensitivity of a directed spin- 
down search would be to constrain the spindown parame- 
ter space by taking advantage of whatever information we 
have about the source from physical models. In practice 
we do not have and probably will never have completely 
reliable models for this type of system. The phenomeno- 
logical model presented here has many physical assump- 
tions that may prove to be inadequate. For instance, in 
the saturation phase it is highly unlikely that the star 
can be represented simply as a uniformly rotating equi- 
librium configuration plus a linear perturbation. When 
the mode reaches saturation, the mean velocity perturba- 
tion is comparable to the rotation rate of the star. Thus 
in the saturation phase, the velocity field may develop 
complicated nonlinear structures (such as the cyclones on 
Jupiter) that produce gravitational radiation involving 
many different multipoles. Long-lived non-axisymmetric 
structures requiring higher-order modes are seen in the 
density perturbations of simulations of star formation 
when there is enough angular momentum [0. Although 
the physics is very different in such simulations, a priori 
we see no reason to believe that the velocity analogues of 
these structures are not formed by the r-mode instability. 

However, not all of the physics affects the waveforms. 
Some kind of phenomenological model of the signal (as 
opposed to the star) could be enough to substantially 
reduce the volume of spindown parameter space to be 
searched. For a reasonable model of the neutron star, 
all of the terms in the expansion for df /dt might be de- 
termined from a relatively small number of phenomeno- 
logical parameters. Although quite crude, the model of 
r-mode instability gravitational waveforms provided here 
has several features which should be fairly robust: the 
form of h(f) during the spindown phase, the approximate 
frequency range of the expected radiation, the approxi- 
mate timescale for the spindown to occur, etc. Presum- 
ably these robust features can be used to reduce consid- 
erably the volume of general spindown parameter space 
which need be searched. 
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FIG. 10. Spindown age for an r-mode-driven young neu- 
tron star for different values of the nonlinearity parameter k. 



B. Stochastic Background 

We now consider the gravitational-wave stochastic 
background generated by spin-down radiation from neu- 
tron star formation throughout the universe. A stochas- 
tic background is detected by looking for correlations in 
the response of two or more detectors. The sensitiv- 
ity of the network to the background drops rapidly for 
gravitational wave frequencies much higher than the in- 
verse light-travel time between detectors. For the present 
application, the important networks are therefore the 
VIRGO-GEO pair (high frequency cut-off / cut w 400 Hz) 
and the Washington-Louisiana pair of LIGO detectors 
(/cut ~ 100 Hz). The gravitational radiation generated 
by the r-mode spin-down process has significant power 
at these relatively low-frequencies, / rs 100 — 400 Hz. In 
addition the cosmological redshift (with z 1 — 4) of 
these sources has the beneficial effect of shifting much of 
the radiation into the detectable band. 

Neutron stars have presumably been formed since the 
beginning of star formation. If each neutron star for- 
mation event converts a reasonable fraction of a solar 
mass into gravitational radiation via the r-mode instabil- 
ity, then the sum of this radiation constitutes a random 
background that may be detectable by LIGO III. We now 
make a rough estimate of the spectrum and detectability 
of this background radiation. A more detailed analysis is 
being carried out by Vecchio and Cutler |Q . 

The spectrum of the gravitational wave background 
is typically represented by the following dimensionlcss 
quantity: 



1 dp gw 
Pc d log / ' 



(5.16) 



where p gw is the energy density in gravitational waves, 
and p c = 3c 2 H$/3nG w 1.6 x 10~ 8 /i 2 00 erg/cm 3 is the 
critical energy density just needed to close the universe. 
(H is the Hubble constant and hioo is H divided by 100 
km/s/Mpc.) The signal-to-noise with which this back- 



ground can be detected in a correlation experiment be- 
tween two detectors (here assumed to have uncorrelated 
noise) is given by [ p3]]34| ] 
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where T is the integration time, 1 Sf l (f) and 2 Sh(f) are 
the noise spectral densities of the two detectors, and 7(/) 
is the (dimensionless) overlap reduction function, which 
accounts for the fact that the detectors will typically have 
different locations and orientations. 

We can get a rough estimate of J7 gw (/) due to the 
stochastic background of gravitational radiation from the 
r-mode instability as follows. A 1.4M Q neutron star ro- 
tating with Kcplerian angular velocity Qk 
rotational kinetic energy Ek ~ 0.025M Q c 
which is radiated as gravitational waves. For now, as- 
sume that all neutron stars are born with angular velocity 
£Ik- For a single supernova occurring at z = 0, the spin- 
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down radiation has spectrum dE gw 
where / max « 2Q K /3n w 1400 Hz. 

We are chiefly interested in 17 gw (/) for / < / min 120 
Hz, since that is the range where the pair of LIGO II 
detectors will have their best sensitivity. Let n(z) dV dz 
be the number of supernovae occurring within co- moving 
volume dV and redshift interval dz. In the frequency 
range of interest then the spectrum of gravitational radi- 
ation in the universe today is given by 
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(5.18) 



where /'A/' = (l + z) 2 /A/ are the values of the frequen- 
cies as emitted by the source, z m - m (f) = max{0, /min// - 
1}, w(/) = nun{z*,/ max // - 1}, and z* corresponds 
to the maximum redshift where there was significant 
star for mati on. The factor 1 + z in the denominator 
in Eq. ( 5.1 8| ) accounts for the redshift in energy of the 
gravitational radiation. 

To evaluate the integral in Eq. ( |5.18 ), we must make 
some assumption about the rate of Types lb, Ic, and II 
supernovae (which are the ones that leave behind neutron 
stars). The combined rate in our galaxy at present is 
roughly one per 100 years. If this rate had been constant, 
the Galaxy would today contain about 10 s neutron stars. 
However at earlier times, for z between 1 and 4, the rate 
(with respect to proper time) was significantly higher. 
A reasonable estimate is that our Galaxy contains 3 x 
10 8 neutron stars today. Let n(z)Az be the density of 
neutron star births (per unit comoving volume) between 
redshifts z and z + Az. As a rough, first-cut at this 
problem, we model n(z) as constant n(z) = n Q for < 
z < z*, where w 4, and ta ke n( z) — for z > z*. In 
this case the integral in Eq. ( 5.18| ) can be performed to 
obtain: 
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^gw (/) 



0, 

(z* + I) 2 - (/„ 

z*(z* + 2), 

(/max/ /) — lj 

o, 



>//r 



< / < fx, 

fl < f < /min) 
/min < f < f%> 
fi<f< /max, 
/ 5* /max. 



where A = 2n o ^ K /(3y0 c / max 

/2 = /max/(l + 2*)- 

We can estimate the value of the constant n Q as fol- 
lows. We assume that the number of neutron stars in 
a given location is roughly proportional to the lumi- 
nosity of the visible matter at that location. The to- 
tal luminosity of our Galaxy is 1.4 x lO 10 ^© |3^], while 
the number of neutron stars in the Galaxy is about 
3 x 10 8 . Thus the neutron star mass to total luminos- 
ity of matter ratio is about Q.03M & / Lq. The mean lu- 
minosity of the universe is 1.0 x 10 8 /ii o-k©/Mpc 3 [j36| . 
Thus, the mean mass density of neutron stars is about 
p ns w 3 x lO 6 /ii OO M /Mpc 3 w 1.1 x 10- 5 hio p c . The 
current density of neutron stars is related to the con- 
stant n by p ns « 1.4M z*7i o w 5.6M n o . Thus 
^n E K /p c ~ 3.3 x I0~ s h^n- For LIGO, the most impor- 
tant range in Eq. ( |5. 19 ) is f\ < / < / m i n ; we can re- write 
the result for this portion in the more useful form 



(5.19) 

), /i = /min/(l + z*) and 



Ogw(/) » 2.4 x \Q- W h 



100 



/ 



24 Hz 



- 1 



(5.20) 



Evaluating at / = 50 Hz, with h\oo = 0.7 we find 
fj gw (/ = 50 Hz) w 1.1 x 10~ 9 . 



Using the above spect rum Eq. ( |5.19| ), we evaluate 
the S/N using Eq. ( f5.17| ). For an integration time of 
T = 10 7 s we find SjN = 0.0022, 0.34, and 2.6 for 
LIGO I, II, and III respectively. Since there has been 
some discussion of building a second kilometer-size in- 
terferometer in Europe, we also consider the sensitivity 
to the r-mode background of this detector paired with 
VIRGO. We assume that the detectors will be located 
less than about 300 km apart and have the same orienta- 
tion. (To model this, we simply set j(f) = 1 in the inte- 
gral Eq. [5.17 ). We find S/N = 5.6 assuming both these 
detectors have LIGO II sensitivity, S/N = 0.9 assuming 
one detector has LIGO I sensitivity while the other has 
LIGO III sensitivity, and S/N — 23 assuming both have 
LIGO III sensitivity. Thus we see that detection of the 
r-mode background will have to wait either for develop- 
ment of "advanced" interferometers or for the construc- 
tion of two nearby detectors with "enhanced" sensitiv- 
ity. Two nearby "advanced" interferometers could see 
quite a strong signal. All the above results on correla- 
tion measurements of the stochastic background assume 
that magnets will be eliminated from the LIGO design. 
With the current design long-range correlated magnetic 
fields from Schuman resonances and lightning strikes will 
mimic a stochastic background with J7 gw approximately 

10~ 7 to io~ 9 §71. 



These calculations assumed that all neutron stars are 
born with spins near their maximal value fix ■ It should 
be clear, however, that these results for the S/N achiev- 
able by the LIGO pair depend only on the stars being 
born with spins greater than about 300 Hz. Of course, it 
could well be that some fraction F of neutron stars are 
born with rapid spins, while (1 — F) are born slowly spin- 
ning. The values of S/N in this case could be estimated 
from those give n above by multiplying by F. (See Spruit 
& Phinney ]38[ for a recent argument that most neutron 
stars should be born with very slow rotation rates.) 

It has previously been suggested that there could be 
a detectable gravitational wave background produced by 
supernova events |39 40 1 . The stochastic background due 
gravitational radiation from the r-modes differs from that 
previously envisioned in two important respects. Previ- 
ously it was assumed that the radiation would be emit- 
ted in short bursts, forming a random but not continu- 
ous background. For the r-mode background, the long 
duration of the emission guarantees that it will be a con- 
tinuous hum rather than an occasional pop. Also, the 
spectrum from spindown radiation extends to lower fre- 
quencies than had previously been expected from super- 
nova events. 



VI. DISCUSSION 

The discovery of a strong source of gravitational waves 
that is ubiquitous and is associated with such interesting 
objects as supernovae and neutron stars inevitably opens 
up a rich prospect for obtaining astronomical information 
from gravitational wave observations. We shall discuss 
here some of the more obvious prospects. 

Background radiation from r-modes. Pleasantly, the 
background requires no detailed modeling of the signal 
in order to detect it. However, detection of background 
radiation from r-modes will probably have to wait for 
the development of "advanced" interferometers (or the 
construction of two nearby "enhanced" interferometers), 
even if we assume that a large fraction of neutron stars 
are born rapidly rotating. For nearby detectors with 
LIGO III sensitivity, S/N is high enough that one could 
experimentally measure the spectrum £l gw (f) with rea- 
sonable accuracy. This might provide very interesting 
cosmological information. For instance, imagine that 
most neutron stars are born rapidly rotating. (This could 
be verified, at least for neutron stars born today, by di- 
rect LIGO detections of nearby supernovae.) The back- 
ground spectrum between 25 and 50Hz would then give 
us direct information about the star formation rate in 
the early universe. And its spectrum between 50 and 
perhaps 300 Hz would tell us about the distribution of 
initial rotation speeds of neutron stars. 

Individual r-mode events associated with known super- 
novae. Observations of individual spindown events will 
be more difficult to achieve but can be very rewarding. 
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The easiest case is if the supernova that leads to the neu- 
tron star is seen optically. This gives some hint of when 
the r-mode radiation should be looked for, but more im- 
portantly it gives a position. That reduces the difficulty 
of extracting the signal from the detector data stream. 
Detection of the radiation will return the amplitude of 
the signal, its polarization, the final spin of the star, and 
the values of the parameters of the waveform. Assuming 
that the three large detectors (the two LIGO installations 
and VIRGO) all detect the signal with S/N « 8, the ef- 
fective combined S/N will be 8y/3 ss 14. Values of the 
various parameters will then typically be determined to 
10-20% accuracy. The polarization measures the orien- 
tation of the spin axis of the neutron star, which will be 
difficult to relate to any other observable. So it may not 
be of much interest. But the power spectrum of the radi- 
ation measures, as we have seen, the loss of rotational en- 
ergy of the star, so it is proportional to I/D 2 . If the host 
galaxy's distance can be determined to better than the 
accuracy of the gravitational wave measurement, which 
seems likely, then this will provide a direct measure of 
the moment of inertia of a neutron star. 

There will likely be several detections per year, which 
will shed light on a number of uncertainties. Even if the 
parameters are only the Taylor expansion coefficients for 
the frequency, they will constrain models of the r-mode 
spindown. We can expect to get some information about 
cooling rates, viscosity, crust formation, the equation of 
state of neutron matter, and the onset of superfluidity 
(or some combination of these). We also expect variabil- 
ity from event to event, due to different initial conditions 
after gravitational collapse, such as differential rotation 
or even the mass of the neutron star. Significant differ- 
ential rotation might affect the final spin rate of neutron 
stars; hence any variability in the final spin speed might 
shed light on these initial conditions. 

If we find we can detect this radiation with confidence, 
then the absence of it after a supernova could be a hint 
that a black hole has formed instead of a neutron star. 
More rarely it might indicate that the neutron star re- 
mained bound in a binary system, and the orbital mod- 
ulation of the signal made it impossible to find. 

Some supernovae will be especially nearby, and some 
neutron stars will be oriented more favorably, so there 
might be a handful of events with single-detector S/N w 
15 — 20. These will provide particularly good constraints 
on the moment of inertia, superfluidity, viscosity, etc. For 
these events it may also be possible to trace the waveform 
back to its initial stages, and thereby to measure the 
initial spin of the collapsed star. In some cases, a star 
might be formed with more mass that it can support after 
spinning down, and so the r-mode spindown would lead 
to and be interrupted by collapse to a black hole. This 
would probably happen relatively soon after the star is 
formed; but for strong events, it might be possible to 
detect this break. 

Using r-mode events as supernova detectors. Perhaps 



the most exciting use of r-mode observations would be to 
identify hidden or unnoticed supernovae. If it proves pos- 
sible to create search strategies that are efficient enough 
to detect r-mode spindown even without prior positional 
information from an optical observation, then the gravi- 
tational wave detectors will become supernova monitors 
for the Virgo cluster. Perhaps as many as half of the 
supernovae in Virgo go unnoticed, hidden in thick dust 
clouds. LIGO and VIRGO would not give optical ob- 
servers advance notice of the supernova: they will iden- 
tify a neutron star only a year or so after it was formed. 
But they may be able to locate the position of the event 
with an accuracy of better than one arc-second. 

This great precision is achieved from the modulation 
of the signal produced by the motion of the earth |fll| . 
The angular accuracy is similar to that achieved for pul- 
sar observations in the radio. Fundamentally it is the 
diffraction limit of a gravitational wave telescope with 
the diameter of the earth's orbit, because the detector 
acts like a synthesis array as it builds up signal along 
its orbit around the Sun. The ratio of the gravitational 
wavelength of about 1000 km to the diameter of 2 AU is 
about 3 x 10 -6 radians. This angular accuracy improves 
with S/N as well, but it can be degraded by uncertainties 
in the spindown parameters. This assumes, as we have 
in this paper, that neutron-star cooling takes a year or 
more. If the alternative cooling scenarios are correct and 
the star cools in a few days, then the angular accuracy 
of observations will be very poor. 

For an event at 20 Mpc, arc-second angular resolution 
corresponds to distance resolution of about 60 pc. Obser- 
vations would therefore not only tell us in which galaxy 
the event occurred, but even in which molecular cloud. 
Detailed follow-up searches will then be possible for the 
expanding nebula, starting perhaps one year after the 
event. 

Other implications for gravitational radiation from su- 
pernovae. The r-mode instability also has implications 
for our expectations of other kinds of gravitational radi- 
ation from supernova events. Although detecting super- 
novae was the goal of initial bar detector development, 
it has not been possible before now to provide reliable 
predictions of radiation from gravitational collapse. The 
r-mode instability is a reliable prediction, but only of 
radiation long after the collapse event. It is clear that 
a collapse that produces a rapidly rotating star will be 
more likely to radiate strongly, especially if it can reach 
the dynamical bar-mode instability that is seen in the 
lower-density star-formation simulations. But although 
rapid rotation is in some sense natural in gravitational 
collapse, the fact that young pulsars like the Crab are 
slow rotators seemed to indicate that neutron stars do 
not form with fast spins. Now the r-mode instability has 
provided an explanation for the slow spins of young pul- 
sars; there is no longer any observational restriction on 
the initial spins of neutron stars. 

It therefore seems to us much more likely than be- 
fore that the gravitational collapse event can also be a 
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strong source of gravitational waves. The following sce- 
nario seems plausible in at least the extreme cases where 
rotation completely dominates the last stages of collapse. 
The collapsed object has so much spin that it forms a bar 
shape on a dynamical timescale. This radiates away an- 
gular momentum in gravitational waves until the star is 
finally able to adopt a stable axisymmetric shape. The 
strong gravitational radiation ceases, to be replaced by 
the developing r-mode radiation. The first burst would 
be detectable by LIGO II at the distance of the Virgo 
cluster, and a network of such detectors could give a 
rough position, which would allow notification of optical 
astronomers of the event and multi-wavelength follow-up 
observations. In the gravitational wave data stream, in- 
tensive searching for the r-mode radiation would follow. 

The unexpected strength of the r-mode gravitational 
wave instability in young neutron stars has therefore 
completely changed the prospects for detection of gravi- 
tational waves from supernova events. Not only can we 
look forward to regular detections of the spindown radi- 
ation when enhanced detectors begin operating, but we 
also now have more reason to expect strong bursts from 
the collapse events themselves. Particularly exciting is 
the prospect for using networks of LIGO-type detectors 
as supernova monitors, able to pinpoint the positions of 
hidden supernovae with arc-second precision. In order 
to achieve these goals, however, much work needs yet 
to be done to eliminate the uncertainties in our models 
of young neutron stars and to develop effective ways of 
searching for signals of this kind. 
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